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Abstract 

During their asymptotic giant branch, evolution low-mass stars lose a significant fraction of their mass through an intense wind, enriching the 
interstellar medium with products of nucleosynthesis. We observed the nearby oxygen-rich asymptotic giant branch star IK Tau using the high- 
resolution HIFI spectrometer onboard Herschel. We report on the first detection of H, 6 and the rarer isotopologues H^O and Hj s O in both the 
ortho and para states. We deduce a total water content (relative to molecular hydrogen) of 6.6 x 10~ 5 , and an ortho-to-para ratio of 3:1. These 
results are consistent with the formation of H 2 in thermodynamical chemical equilibrium at photospheric temperatures, and does not require 
pulsationally induced non-equilibrium chemistry, vaporization of icy bodies or grain surface reactions. High-excitation lines of 12 CO, CO, 
28 SiO, 29 SiO, 30 SiO, HCN, and SO have also been detected. From the observed line widths, the acceleration region in the inner wind zone can be 
characterized, and we show that the wind acceleration is slower than hitherto anticipated. 

Key words. Line: profiles, Radiative transfer, Stars: AGB and post-AGB, (Stars): circumstellar matter, Stars: mass loss, Stars: individual: IK Tau 
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1. Introduction 

IK Tau, also known as NML Tau, is an extremely red, oxygen- 
rich, Mira-type variable, with a period of about 470 days 
dWing & Lockwoodll973l) . Its dust-driven wind produces a cool 
circumstellar e nvelope (CSE), w hich fosters a rich gas-phase 
chemistry (e.g., iDuari et all 1999I). IK Tau is relatively nearby, 
at a distance of ~265 pc dHale et alj |1997 | ). Estimates of its 
mass-loss rates range from 3.8 x 10~ 6 (|Neri et al.l Il998l) to 
3 x 10~ 5 M o /yr (iGonzalez Delgado et al.ll2003l) . IK Tau's prox- 
imity and relatively high mass-loss rate facilitate the observa- 
tion of molecular emission lines. Currently, a dozen different 
molecules and some of their isotopologs have been discovered 
in I K Tau, including CO, HCN, SiO, SiS , SO, S0 2 , and NaCl 
(e.g. lMilam et al.ll2007tlDecin et al.ll2010h . 

In this Letter, we report on the detection of thermal emis- 
sion of water (H2O) in the envelope around IK Tau. The main 
isotopolog (H^O) as well as the rarer isotopologs (H^O and 
H' s O) are detected for both the ortho- and para-states. We also 
present observations of high-excitation rotational transitions of 
12 CO, I3 CO, 28 SiO, 29 SiO, 30 SiO, HCN, and SO and demon- 
strate that the observed line widths characterize the acceleration 
region in the inner wind zone. 



2. Observations and data reduction 



observe molecular fingerprints in the frequency ranges of 480- 
1150GHz and 1410-1910GHz at a spectral resolution up to 
125 kHz. Single-point observations towards IK Tau were carried 
out with the Herschel/HIFI instrument in the dual beam switch 
(DBS) mode with a 3' chop throw. The observation strategy and 
data-reduction are disc ussed in the online Appendix A and in 
lBuiarrabaleTai1d2010l) . 



3. Results 

In the 13 single-point observations obtained so far, covering a 
frequency range of 92.8 GHz in total, 31 molecular emission 
lines have been detected belonging to 12 CO, 13 CO, H 16 0, H' 7 0, 



The HIFI instru ment (Ide Graauw et all 1201 Oh onboard the 
Herschel satellite dPilbratt et al.l 120101) offers the possibility to 



* Herschel is an ES A space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 



Hf O, 28 SiO, 2y SiO, JU SiO, HCN, SO, and NH 3 (see Table |A~21 
in the online Appen dix). The detection of NH3 is described in 
lMentenetai1d2010l) . 

For the first time, different excitation lines of water are dis- 
covered for the two nuclear spin isomers, ortho- and para-H lfi O, 
as well as for the rare isotopologs H^O and H^O (see Fig.[TJ. 
The observations of these transitions provide information about 
the total water content, the ortho-to-para ratio, and the isotopic 
ratios H* 6 0/H 17 and H 16 0/H' 8 (see Sect. ED. 

High-excitation rotational transitions are observed for dif- 
ferent molecules (see Fig. |2). From the observed line widths, 
which range between 11 and 19km/s, it is immediately clear 
that the HIFI observations offer us a strong diagnostic to trace 
the wind acceleration zone in the inner envelope (see Sect. 13.2b . 
Moreover, the high-excitation 12 CO J=10-9 and J=16-15 lines, 
complemented with ground-based observations of CO for 
J=l-0 to J=7-6, can be used as temperature indicator for the 
envelope as close as ~20R+ (see Sect. 13.21 ). 
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Figure 1. HIFI observations of water lines (gray) compared to 
their line profile predictions (black). Row 1-2: ortho-H} 6 0, row 
3-4: para-H^O, row 5: ortho and para-H^O, and row 6: ortho 
and para- Hi, 8 



3.1 . Water in IK Tau's envelope 

Water content: We observed both ortho and para-water lines 
for different isotop ologs in IK Tau ( Fig.[T|). Using the analytical 
formula from iGroenewegenl dl994l) . we calculated a photodis- 
sociation radiu s of 2.8 x 10 16 cm or 187 0/?*, while the theoret- 
ical models by IWillacv & MillaJ d 1997b predict a value around 
1600/?+. Assuming the same photodissociation radius for all iso- 
topologs and isomers, we deduced (see Sect.|B]in the Appendix) 
a water-abundance [ortho-H2 6 0/H2]=5 x 10~ 5 and an ortho-to- 
para water ratio (OPR) of 3: 1 . Taking the uncertainty in the pho- 
todissociation radius into account, the estimated uncertainty in 
the water abundance is a factor 2 for the assumed envelope struc- 
ture (see Sect. 13. 2t . 

The total (ortho+para) H^O abundance for IK Tau derived 
in this Letter, relative to the total hydrogen content and assum- 
ing that all hydrogen is in the form of molecular hydrogen, is 
[H 2 6 0/H tot ]=3.3 x 10" 5 , slightly lower than the photospheric 
thermodynamic equilibrium (TE) predi ction of 7 x 10~ 5 for an 
evolved star with a C/O ratio of 0.75 dCherchnefHl2006l) . The 
theoretical TE value for the total water abundance is, however, 
very dependent on t he exact value of the C/O ratio, which is 
unknown for IK Tau. ICherchnefn d2006l) predicts a higher inner 
wind abundance of 3.5 x 10~ 4 around 5 R+ in the case of pul- 
sationally induced non-equilibrium chemistry. The upper limit 
implied by the cosmic abundances of carbon and oxygen is no 
more than 1 x 10~ 3 , assuming that all carbon is locked in CO an d 
the remaining oxygen goes i ntoH?0 dAnders & Grevessell989t) . 
Vaporization of icy bodies dJusttanont et alj|2005l) or grain sur- 
face reactions via Fischer- Tropsch catalysis dWillacvll2004l) will 
also result in a higher water abundance. Considering the agree- 
ment between the derived and the TE abundance, the effects of 
the other processes seem negligible for IK Tau. 



On the basis of ISO-LWS data, iMaercker et"aD d2008l) de- 
rived for IK Tau an ortho-H^O abundance (relative to H2) of 
3.5 x 10~ 4 , significantly higher than the value we derived. Both 
studies used an almost equal photodissocation radius. However, 
IMaercker et al.l d2008l) only inclu ded the ground-sta te for the 
excitation analysis. As shown by iDecin et all d2010t) . neglect- 
ing the first vibrational state of the bending mode, vn — 1, 
leads to an underprediction of the HIFI transi tions here pre- 
sented. Since the ISO data used in the study of IMaercker et al.l 
are saturated, the integrated intensity is insensitive to the exact 
value of the abundance, which will result in an overprediction 
of the H20-abundance. We also note that the H2O line pro- 
file predictions are quite sensitive to the exact value of the ki- 
netic temperature and the dust radi ation field (and w avelengfh- 
dependent absorption efficiencies) dDecin et al.ll2010l) . which in 
part might also explain the difference between t he value we de- 
rive [ortho-H' 6 0/H 2 ]=5 x 10~ 5 and the value of IMaercker et al.l 
d2008h [ortho-H 2 6 0/H 2 ]=3.5x 10~ 4 . A comparison with the ISO- 
LWS data will be presented in Decin et al. {in prep.). 

For the oxygen-rich Mira W Hva. IZubko & Elitzuil d2000h 
derived an OPR value of 1:1.3 and Barlow et al.l d 19961) reported 
a value of 1:1. However, iBarlow et alJ noted that their derived 
value is q uite uncertain, due to th e high opacities in the water 
lines used. Justtan ont et alJ (l2010h derived an OPR of 2. 1 : 1 for 
the S-type AGB x Cyg. Our study here is the first time that the 
OPR value in an oxygen-rich Mira is determined from a com- 
bination of emission lines of H' 6 0, H^O, and H' s O. The ad- 
vantage of using the rarer isotopologs is that the opacity of the 
lines is lower, making the line intensities more sensitive to the 
OPR value. The observed para-water lines are consistent with 
an OPR value of 3:1 (+0.4). The lowest energy level of para- 
H2O is ~34 K below that of ortho-H20 When water forms in the 
gas phase via exothermic reactions the energy released is much 
greater than this energy difference and the OPR reflects the high- 
temperature (~50K) thermodynamic 3:1 ratio of the statistical 
weights between the species. The derived OPR value of 3: 1 con- 
firms that water in IK Tau is formed in warm and dense regions 
of the envelope where the chemistry is in thermodynamical equi- 
librium. 



Isotope ratios: Assuming the same photodissociation radius as 
for the main isotopolog, the isotopic ratios we derive for IK Tau 
are H} 6 O/H' 7 O=600 (±150) and H 2 6 O/H} 8 O=200(+50), hence 
well below the sola r values ( 16 0/ 17 0~2632 and 16 0/ 18 0~499; 
lAsplund et al.l 120091) . Interpreting the derived isotopic ratios in 
terms of nucleosynthesis and subseq uent dredge-ups or extra 
mixing processes is quite complex (e.g. lHarris et al.fl 9851, 1 19871 : 
iKarakas et al.ll2010l) . In stars that are sufficiently massive to un- 
dergo CNO-cycle hydrogen burning, the low initial n O abun- 
dance (assumed to be solar) is enhanced. When helium burn- 
ing begins inside the hydrogen-burning shell, 17 is expected 
to be completely destroyed in the region where maximum hy- 
drogen burning occurs. The isotope ls O, on the other hand, is 
expected to be destroyed during hydrogen burning, so that it vir- 
tually disappears from the hydrogen-burning zone and from the 
hydrogen-exhausted CNO equilibrium zone within it. When he- 
lium burning starts, the 18 abundance might slightly increase. A 
succession of convective mixings brings to the surface material 
that is affected by t hese nuclear transformations. Calculations by 
iHarris et all d 19851) show that in every star that becomes a red gi- 
ant star (M£0.8M Q ), the initial 16 0/ 17 ratio decreases to -440 
during the first dredge-up, while the 16 0/ 18 slightly increases. 
The second dredge-up occurs at the end of core helium burning 
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only for the most massive intermediate-mass stars (M^4.5 M ). 
The estimated 16 0/ I7 ratio ranges between 150 and 500, while 
the 16 0/ 18 ratio slightly increases. The third dredge-up occurs 
in the subsequent helium shell-burning phase for stars >2M S , 
and is expected to yield isotopic ratios of 16 0/ 17 < 200. If hot 
bottom burning occurs (for stars above 3-4 M Q ), the 16 0/ 17 ra- 
tio will be of the order of 20-50. That the 16 0/ 17 ratio is around 
600, implies that the first but no subsequent dredge-ups occurred 
and constrains the initial mass of IK Tau to be within 1-2 M . 
Alternatively, if the star is more massive than 2 M Q and the third 
dredge-up has occurred (but has not turned the star into a carbon- 
rich star), transferred material from a post-third dredge-up enve- 
lope must have had a 16 0/ 17 ratio < 200. This implies that in 
stars with 16 0/ 17 ~ 600 the transferred material has been heav- 
ily diluted by material from the star's own envelope with much 
higher 16 0/ 17 ratios. From Fig. 5 of lHarris et ail dl987l) . it is 
estimated that IK Tau has a low s-process neutron exposure, To 
< 0. 1 , implying a low absolute enhancement of the s-process el- 
ements and only a few third dredge-up events, consistent with 
IK Tau still being an oxygen-rich AGB star. 

Hitherto, the lower than solar 16 0/ 18 ratios cannot be ex- 
plained by any stellar evolution model in the literature. However, 
IK Tau is not the only Galactic star with a low 16 0/ 18 ratio (see 



Fig. 3 in Karakas et al.ll20Tob ; some barium stars analyzed by 



lHarris et all (1 19851) also have a low 16 0/ 18 value. It is antici- 
pated that the observations of other evolved stars in the frame- 
work of the HIFISTARS programme (P.I. V. Bujarrabal) will add 
new information to this discussion. 

3.2. Thermodynamics! structure of the envelope 

Being a simple diatomic molecule with a well understood en- 
ergy diagram, CO has been successfully use d to study the struc- 
ture of the CSEs around evolved stars (e.g. JSchoier et al1 l2002: 
iDecin et ai]|2006l) . Different transitions can be used to investi- 
gate different regions of the envelope, prob ing the dens i ty, the 
temperature, and the velocity of the CSE. IDecin et aU J2010I) 
used the 12 CO J=l-0 to J=7-6 lines to determine the thermo- 
physical structure of the CSE of IK Tau beyond ~100/?*, based 
on a non-local thermodynamic equilibrium (non-LTE) radiative 
transfer analysis of the available transitions. In the first instance, 
the kinetic temperature and velocity structure of the envelope 
were calculated by solving the equations of motion of gas and 
dust and the energy balance simultaneously. To get insight into 
the structure in the inner wind region, the HCN J=3-2 and J =4- 
3 transitions were used, since observational evidence exis ts that 
HCN is formed close to the star (<3.85", iMarvell 120051) . The 
Gaussian HCN line profiles indeed point toward line formation 
partially in the inner wind where the stellar wind has not ye t 
reached its fu ll terminal velocity dBuiarrabal & Alcolealll99ll) . 
The results of IDecin et alj d2010l) infer a wind acceleration that 
is lower than derived from solving the momentum equation. 

3.2.1. Temperature structure 



Adopt ing the thermodynamic structure derived in IDecin et ail 
d2010h (and reproduced in Fig. IB.ll of the online Appendix), the 
theoretical line profiles for the 12 CO J=10— 9 and J=16-15 are 
calculated (see Fi g. [2]l using the non-LTE ra diative transfer code 
GASTRoNOoM dDecin et alj I2006L l2010h . The 12 CO J=10-9 
line is very well reproduced, while the J=16-15 line exhibits 
slightly larger deviations. The latter most likely reflects the very 
difficult calibration of this frequency setting (for which stand- 



ing w aves heavily perturbed the baseline, see iBuiarrabal et al.l 
2010). This result is consistent with the temperature stru cture i n 
the region between 20 and 100/?* derived bv lDecin et alJd2010t) . 

3.2.2. Velocity structure 

To constrain the wind acceleration in the CSE, all molecular 
emission lines as shown in Figs. QJ[2] were modelled (see the 
online AppendixIBl. The line formation region of each molecu- 
lar line was estimated by considering the range of projected radii 
where I Vo (p) pdp, with I Vo the intensity at the line center and p 
the impact parameter, exceeds half of its maximum value (see 
Fig. [3]). We note that the radial extent of the line formation re- 
gion is almost insensitive to the exact value of the velocity struc- 
ture. Optical depths effects can strongly affect the observed line 
widths, and detailed radiative transfer modelling (as presented 
in this Letter) is required to determine the underlying velocity 
structure. Most of the observed lines have line widths in excess 
of 17 km/s. However, a few lines are considerably narrower, and 
their line formation regions are located in that part of the en- 
velope where the wind has not yet reached its terminal velocity 
(19 km/s). Although the line formation regions are quite broad, 
we find the first observational evidence that the wind acceler- 
ation is slower than im plied by the mome ntum equation, cor- 
roborating the results of Decin et all d2010t). Using the classical 
/?-parametrization (e.g.. lLamers & Cassinellil[l999t) to simulate 
the velocity structure, 



v(r) - u + (uoo - u )(l -I 



(1) 



where uq is the velocity at the dust condensation radius and v^ is 
the terminal velocity, we find that 1 < ji < 2. The theoretical line 
predictions shown in this Letter are computed for/? = 1 . The mo- 
mentum equation, in contrast, matches a much steeper velocity 
profile with /?=0 . 6. We note that the velocity structure derived by 
lJusttanont et al.l d2010|) for y Cyg is com pliant with a/?- value of 
0.9. As discussed bv lDecin et alj d2010l) . there are several possi- 
ble causes for a less steep velocity structure. We summarize that 
a slower wind velocity may be caused by incomplete momentum 
coupling, dust emission being slightly optically thick to the stel- 
lar radiation, the fact that not all dust species are formed at the 
same time and at the same radial distance, and/or that some dust 
species are inefficient as wind drivers. 



4. Conclusions 

Using the HIFI spectrometer onboard Herschel, we have ob- 
served the evolved oxygen-rich Mira star, IK Tau, in 3 1 molec- 
ular emission lines. For the first time, several lines of ortho and 
para H' 6 0, and the rarer isotoplogs Hj, 7 and Hj 8 have been 
detected. We have deduced a total water content (relative to H2) 
of 6.6 x 10~ 5 , and an ortho-to-para ratio of 3:1. The measured 
thermal emission from water therefore points toward TE chem- 
istry, as opposed to pulsationally induced non-equilibrium chem- 
istry, grain surface reactions, or evaporation of icy bodies. For 
the isotopic ratios, we find subsolar values H' 6 O/H' 7 O=600 and 
H' 6 O/H' 8 O=200. 

The high-excitation lines of 12 CO, 13 CO, 28 SiO, 29 SiO, 
30 SiO, HCN, and SO indicate that the wind acceleration is slower 
than hitherto anticipated from solving the momentum equation. 
These results show the great capability of HIFI to study the com- 
plex thermodynamical and chemical envelopes of evolved stars 
in tremendous detail. 
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Figure 2. HIFI observations of line transitions of several molecules (gray) compared to their lin e profile prediction s (black). Full 
black lines refer to line predictions using the molecular abundance stratifications as determined bv lDecin et alJd2010l) . dotted black 
lines show the model predictions for the refined fractional abundances as determined in this Letter (see the online Appendix B). 
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Figure 3. Velocity profile of IK Tau. Velocity data are ob tained 
from mappin g of maser emis sion: SiO dBoboltz & Diamon 






l200l . H 2 (lBainsetal.H200l . and OH ([Bowers et al.l ll989 
The CO expansion veloc ity derived from gr ound-based CO J=l- 
data is also indicated dDecin et al.ll2010h . The triangles show 
the place in the envelope where the line formation is highest for 
the HIFI data presented in Figs. QJ[2] i.e., where I Vo (p) p 2 is at 
its maximum. The horizontal bars show the minimum and max- 
imum radial distance for the line formation of each individual 
transition. The vertical bars show the uncertainty on the observed 
line widths. The expansion velocity deduced from solving the 
momentum equation is shown by the full green line. The dashed 
blue line represents a power law (Eq. [TJ with /? = 1 (as used 
for the modelling in the online Appendix iBl. For comparison, 
an even smoother expansion velocity structure with /3- 1.8 is 
shown with the red dotted line. The vertical dashed black line in- 
dicates the dust condensation radius Rinner. The velocity at Ri nne r 
is assumed to be equal to the local sound velocity. 
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Appendix A: Details about the observation strategy 
and data reduction 

The regular DBS mode was used for bands 1 up to 5, while 
FastDBS was applied to bands 6 and 7 to achieve higher sta- 
bility with respect to electrical standing waves. Two orthogo- 
nal polarization s were measured s imultaneously. The double- 
sideband (DSB. lHelmich et al . 2010) observations ensure an in- 
stantaneous 8 GHz and 5.2 GHz frequency coverage by the wide- 
band spectrometer (WBS) for respectively bands 1 up to 5, and 
bands 6 and 7. The spectral resolution is 0.5 MHz. Care was 
taken in choosing local oscillator (LO) frequencies such that no 
strong lines from the two sidebands would blend, and that at the 
same time a maximum number of molecular lines would be cov- 
ered. 

The data were processed with the standard HIFI pipeline us- 
ing HIPE, and non-stitched Level-2 data were exported using the 
HiClass tool available in HIPE. Further processing, i.e. blanking 
spurious signals, baseline removal, stitching of the spectrome- 
ter subbands, and averaging, was performed in CLASS. When 
the quality of the spectra measured in both horizontal and verti- 
cal polarization was good, these were averaged to lower the rms 
noise. This approach is justified since polarisation is not a con- 
cern for the presented molecular-line analysis. In all cases, we 
assumed a side-band gain ratio of one. 

All data presented in this Letter were converted from the 
antenna-temperature (TV) scale to the main-beam temperature 
(Tmb) scale according to 7mb - TI /tjmb, wim values of the 
main-beam efficiency t?mb (Table[AJJi calculated for the LO fre- 
quency vlo, according to 



Table A.2. Observation summary. First column gives the molec- 
ular name, second column the transition, third column the fre- 
quency of the line, fourth column the integrated intensity, fifth 
column the noise, and last columns some comments (if neces- 
sary). One line is still unidentified, and is indicated with a 'U' in 
the first column. 



r]MB = r/B x exp 



Vlo/GHz ^ 2 
6xl0 3 



Molecule 


Transition 


Freq. [GHz] 


I [K kin/s] 


noise [K] 


comments 


o-Hj b O 


1 1,0 - 1 0. ] 


556.936 


10.22 


0.005 




o-h; 6 o 


5 3 ,2 - 44,1 


620.701 


14.89 


0.005 


maser 


ott^O 


3 1,2 - 3 ,3 


1097.365 


14.44 


0.036 




o-h; 6 o 


3i,2-2 2) , 


1153.127 


22.60 


0.062 




o-Hl 6 


3n,3 — 2i,2 


1716.770 


30.89 


0.104 




o-H™0 

o-h; 6 o 


5 3 ,2 - 5 2 ,3 


1867.749 


3.45 


0.118 


standing waves 


32,1 - 3l ,2 


1162.912 


13.70 


0.062 




p-Hi 6 


2i,i - 2 ,2 


752.033 


13.71 


0.015 




P-H? 6 

p-h! 6 o 


5 2 ,4 - 4 3 ,i 


970.315 


19.34 


0.026 


maser 


2o,2 - 1 1,1 


987.927 


27.85 


0.030 




P-H? 6 


1 1,1 - 0o,o 


1113.343 


25.10 


0.028 




p-H^O 


1 1,1 -Oq,o Vi = 1 


1205.798 


0.59 


0.070 


tentative, see Fig. 1 


p-h; 6 o 


42,2 -4,, 3 


1207.639 


8.21 


0.070 


maser 


p-h; 6 o 


53,3 "" 6q,6 


1716.953 


1.03 


0.104 


standing waves 


o-HpO 


3l,2 - 3o,3 


1096.414 


1.20 


0.036 




p-HpO 


1 1 , i - Oo.o 


1107.167 


1.50 


0.036 




o-h; 8 o 

p-Hi?0 


3 1,2 - 3o,3 


1095.627 


2.52 


0.036 




1 1,1 - 0o,o 


1101.699 


3.57 


0.028 




12 C 16 


10-9 


1151.985 


9.61 


0.062 




12 C 16 


16-15 


1841.35 


2.48 


0.03 


standing waves 


13 C 16 


9- 18 


991.329 


2.43 


0.030 




13 C 16 


10-9 


1101.350 


2.75 


0.028 




28 Si 16 


14-13 


607.599 


2.31 


0.005 




29 si 16 Q 


13-12 


557.179 


0.89 


0.005 




29 si 16 


26-25 


1112.833 


0.74 


0.028 


low S/N, see Fig. 2 


30 si 16Q 


26-25 


1099.708 


0.39 


0.028 


tentative, see Fig. 2 


HCN 


7-6 


620.304 


2.65 


0.005 




HCN 


13-12 


1151.452 


2.22 


0.062 


tentative, see Fig. 2 


SO 


13 w -12 n 


560.178 


0.20 


0.005 




NH, 


lo-O, 


572.498 


4.93 


0.005 




U 




988.255 


1.37 


0.030 


or at 1003.215 MHz 



X77F, 



(A.l) 



with rjB-0.72 and /7f=0.96 being the beam efficiency in the Hz 
frequency limit and the forward efficiency, respectively. The ab- 
solute calibration accuracy ranges from 10% for the lowest fre- 
quency lines up to 30% for the high frequency (>1 THz) lines. 



Table A.l. LO frequencies vlo, main beam efficiencies 77MB 
according to Eq. IA.ll velocity resolution Av (for 0.5 MHz fre- 
quency resolution), and operational days (OD) when observa- 
tions were carried out. 



vlo 


?7mb 


Ay 


OD 


(GHz) 




(km/s) 




564.476 


0.69 


0.27 


292 


614.763 


0.68 


0.24 


292 


758.772 


0.68 


0.20 


297 


975.087 


0.67 


0.15 


293 


995.481 


0.67 


0.15 


293 


1102.757 


0.67 


0.14 


294 


1 106.735 


0.67 


0.14 


294 


1157.495 


0.67 


0.13 


296 


1200.717 


0.66 


0.12 


296 


1713.596 


0.64 


0.09 


298 


1757.417 


0.63 


0.09 


298 


1838.039 


0.63 


0.08 


298 


1864.539 


0.63 


0.08 


298 



Appendix B: Radiative transfer modelling 

The molecular emission lines shown in Figs. Q~|-|2] were mod- 
elled using the non-LTE radiative transfer code GASTR0NO0M 
(iDecin et al.l l2006h . Last updates to the code and a discus- 
sion of t he available line l ists and collisional rates can be 
found in IDecin et all d20 1 Ot) . The fhermodyna mical structure 
(see Fig. IB. II ) determined bv IDecin et al.l d2010h was confirmed 
using the new HIFI observations (Sect. 13.2b . Modelling the CO 
and H2O lines provides insight into to the co oling/heating rates 
by tr ansitions of these molecules (see also IDecin et aD l2006. 
120101) . As can be seen in Fig. IB. 31 H2O transitions provide the 
main cooling agent in the region up to ~2xl0 15 cm, while adia- 
batic cooling takes over for the region beyond ~2xl0 15 cm. 

In the first instance, the m olecular abundance stratifications 
derived bv lDecin et al.l ( 120101) were assumed to model the 13 CO, 
28 SiO, 29 SiO, 30 SiO, HCN, and SO lines (see full black lines in 
Fig. [2] and full lines in Fig. IB.2t . Using the new HIFI observa- 
tions, the 13 CO and 28 SiO abundance fractions were refined in 
the inner envelope (see dashed black lines in Fig. [2] and dashed 
lines in Fig. IB.2t . 

The radiative transfer modelling for water included the 45 
lowest levels of the ground state and first vibrational state (i.e. 
the bending mode vo — 1 at 6.3 fim) for all isotopologs. Level 
energies, frequencies, and Einstein A coefficients we re ext racted 
from the HITRAN water line list (iRothman et al.ll2009t). Th e 
H2O-H2 collisional rates were taken from iFaure et all d2007l) . 
The effect of including excitation to the first excited vibrational 
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Table B.l. (Circum)stellar parameters for IK Tau dDecin et al.l 
120101) . T e ff is the effective stellar temperature, R+ denotes the 
stellar radius, M the gas mass-loss rate, Ri nner the dust condensa- 
tion radius, Doo the terminal velocity of the wind, and u tul -b the tur- 
bulent velocity in the wind. The molecular fractional abundances 
are given relative to H tot =n(H)+2n(H2), and denote the abun- 
dance at the dust condensation radius (see Fig. IB. 2b . Values be- 
tween parentheses denote refined inner wind abundance values 
obtained from the new HIFI observations. The fractional abun- 
dances for all water isotopologs and isomers are based on the 
HIFI data presented in this Letter. 



Teff [K] 


2200 


[CO/H tot ] 


1 x 10~ 4 


R* [10° cm] 


1.5 


12 CO/ 13 CO 


14(7) 


M [M /yr] 


8 x 10- 6 


[ 28 SiO/H tot ] 


8 x 10" 6 (4 x 10~ 6 ) 


distance [pc] 


265 


[ 29 SiO/H tol ] 


3 x 10~ 7 


^Mnnei L^*J 


8.7 


[ 30 SiO/H tot ] 


1 x 10~ 7 


v m [kms~'] 


17.7 


[HCN/H lot ] 


2.2 x 10~ 7 


<V,rb [kms" 1 ] 


1.5 


[SO/H tot ] 


1 x 10~ 6 (1 x 10~ 5 ) 


[ortho-H} 6 0/H tot ] 


2.5 x 10- 5 


OPR 


3:1 


H> 6 0/HfO 


600 


H> 6 0/HfO 


200 



state of the asymmetric stretching mode (V3 — 1) was tested, and 
was found to be negligible dDecin et alJl2010h . 

A good agreement was found for the HCN(7-6) line prov- 
ing that the inner abundance fraction [HCN/H tot ] is ~ 2.2 x 1CT 7 . 
The 13 CO J=9-8 and J=10-9 lines are somewhat un derpredicted 
assum ing a 12 CO/ I3 CO ratio of 14 as obtained by iDecin et alJ 
d2010l) . Decreasing this ratio to 7 yields a better fit, but we point 
out that the line profiles are quite noisy. The 13 CO fractional 
abundance was obt ained assuming the same photodissociation 
radius as for 12 CO dMamon et al.lll988l) . However, if the effect 
of less self-shielding of 13 CO (compared to 12 C O) were more 
important than estimated bv lMamon et alJ dl988l) . the photodis- 
socation radius of 13 CO would be smaller, affecting the low exci- 
tation rotational transitions more than the higher excitation lines 
observed by HIFI. Another possibility might be that the velocity 
structure is steeper than the /3 — 1 power law now assumed in 
the region between ~20 and 150/?*, where these high-excitation 
lines are mainly formed. Since a constant mass-loss rate is as- 
sumed, this would imply a lower density in this region, and 
hence a higher 13 CO abundance fraction to produce the correct 
line intensity. 

Only one higher-excitation 28 SiO line has been observed so 
far. The J=14-13 transition indicates that the inner wind abun- 
dance might be a factor 2 lower than deduced by IDecin et aD 
d2010l) . yielding an inner wind abundance of 4 x 10~ 6 , decreas- 
ing to 2 x 1(T 7 around 1807?*. The isotopolog line of 29 SiO(13- 
12) is very well predicted for an inner abundace of 3 x 10~ 7 ; the 
higher excitation J=26-25 line of both 29 SiO and 30 SiO are con- 
sistent with the HIFI observations. This implies an isotopic ratio 
of 28 SiO/ 29 SiOof 13. 

Us ing the abundance pattern determined by IDecin et al] 
(1201 Oh . the SO(13i4 - 12i 3 ) is quite well predicted. 
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Figure B.l. Thermodynamic structure of the envelope of IK Tau 
derived from the 12 CO J=l-0 to J=7-6 and HCN J=3-2 and 
J=4-3 rotati onal line transition s for the stellar parameters given 
in Table IB. ll ( Decin et al.ll2010l) . The start of the dusty envelope, 
Rinner, is indicated by the dotted line. 
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Figure B.2. Fractional abundance stratifications for 12 CO, 13 CO, 
28 SiO, 29 SiO, HCN, SO, ortho-H' 6 0, para-H^ 6 0, ortho-H' 7 0, 
para-H^O, ortho-H^O, and para-H^O For all molecules (ex- 
cept water) , the f ull line represents the results obtained by 
IDecin et al.l (l2010h . For 13 CO and 28 SiO, the new results based 
on the HIFI data are shown in dotted lines. The fractional abun- 
dances for all water isotopologs and isomers are based on the 
HIFI data presented in this Letter. 
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Figure B.3. Cooling and heating rates in the envelope of IK Tau 
due to different processes (for details, see lDecin et al.ll2006l> . As 
can be seen in the plot, both CO and H2O transitions mainly 
cool the envelope, but in certain restricted ranges can heat the 
envelope. 



